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ABSTRACT 

We present a study of the interaction of small amplitude, unsteady, freestream distur- 
bances with a shock wave induced by a wedge in supersonic flow. These disturbances may 
be acoustic waves, vorticity waves, or entropy waves (or indeed a combination of all three). 
Their interactions then generate behind the shock disturbances of all three classes, an aspect 
that is investigated in some detail, our motivation here being to investigate possible mech- 
anisms for boundary-layer receptivity, caused through the amplification and modification of 
freestream turbulence through the shock-body coupling. Also, the possibility of enhanced 
mixing owing to additional vorticity produced by the shock-body coupling is investigated. 
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1. Introduction. 

Boundary-layer receptivity (i.e the growth of instabilities within the boundary layer, 
triggered by some form of external disturbances) has been observed in supersonic wind tun- 
nel experiments in the presence of sound waves radiated by turbulent tunnel- wall boundary 
layers or by freestream turbulence. In such situations, the acoustic wave mainstream dis- 
turbances are first processed by the shock wave induced by the body in supersonic stream, 
and then interact with the boundary layer to cause instability (if at all). In this paper, 
we study the former phase of how a shock in the presence of the body processes main- 
stream disturbances. The latter phase of these disturbances internalizing in the boundary 
layer will be subject of future study. The shock-body combination is shown to produce 
a wide spectrum of wavelengths from monochromatic waves, i. e a freestream disturbance 
with a single wavelength; in particular, we note the production of shorter wavelength 
disturbances than would be expected from the free shock-disturbance interaction. Then 
these shorter wavelength perturbations resulting from the shock-body combination would 
go on to interact with the boundary layer on the surface of the body, leading to possible 
instability. In the case of turbulence-shock-body interaction, it is noted that the turbu- 
lence itself is composed of a continuous spectrum of wavelengths and the interaction would 
produce even more short-wavelength disturbances than the single wavelength disturbance 
considered here. Such interactions could possibly enhance mixing of a multi-component 
flow owing to the amplification and generation of vorticity behind the shock. Again, it 
is found that the shock-body combination efficiently processes the freestream disturbance 
and generates a vorticity profile behind the shock which are larger than those arising from 
the interaction of the free shock with mainstream disturbances. 

To investigate these phenomena, we consider the model problem of a wedge in a 
uniform supersonic stream. In the presence of the wedge, the supersonic flow abruptly 
changes direction through an oblique shock wave, details of which are determined by the 
Rankine-Hugoniot conditions. This paper is concerned with the response of the overall 
flow field to general, small amplitude disturbances in the freestream that are convected 
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through the shock wave. In order to quantify this response, we consider isolated single 
wavelength disturbances of acoustic type, vorticity type or entropy type. In particular, 
we are concerned with the production of a continuous spectrum of wavelengths produced 
in the flow behind the shock in the presence of the wedge. Moore (1954), Ribner (1953) 
and McKenzie and Westphal (1968) have investigated the processing of single wavelength 
disturbances by freely propagating shocks, and note that through the interaction of the 
shock and the freestream disturbance, single wavelength entropy, acoustic and vorticity 
disturbances are all generally produced behind the shock in response to an isolated acoustic 
disturbance, or to an isolated entropy disturbance or to an isolated vorticity disturbance. 
With the presence of the wedge, the flow behind the shock is no longer comprised of single 
wavelength acoustic, entropy and vorticity disturbances. In order to satisfy the condition 
of zero normal velocity at the wedge surface as well as the condition that the shock remain 
attached to the apex of the wedge, a continuous spectrum of wavelengths for the acoustic, 
entropy and vorticity disturbances exists behind the shock. It is this more complicated 
disturbance pattern that goes on to interact with the boundary layer developing along the 
surface of the wedge, possibly initiating new receptivity mechanisms. 

Limited work appears to have been done on the theoretical/ computational approach 
to the interaction between shock waves and boundary-layer instabilities, although most 
of this work has focused on the wedge problem. A fairly restrictive (and somewhat ad 
hoc) model was presented by Petrov (1984), using inviscid linear stability equations, and 
somewhat heuristic conditions on the shock (which was also assumed to lie at the edge of 
the boundary layer). Cowley and Hall (1988) presented an asymptotic model, applicable 
to three-dimensional viscous modes of instability, with appropriately simplified conditions 
applied on the shock surface (derived from the Rankine-Hugoniot conditions), which was 
taken to lie just outside of the boundary layer. Chang, Malik and Hussaini (1990) con- 
sidered the viscous small disturbance equations with parallel mean flow assumption, using 
the full inviscid Rankine-Hugoniot conditions on the shock. 

Here we take the equations governing the gas flow both ahead of and behind the shock 
to be the Euler equations (neglecting the boundary layer on the wedge in the first order of 
approximation) along with the ideal gas law. The no mass flux condition is imposed on the 
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wedge surface, so that the normal component of the flow must have the same velocity as 
the velocity of the wedge normal to its surface. Furthermore, it is required that the shock 
remains attached to the apex of the wedge. The analysis is accomplished by linearizing 
the Euler equations about the base state and applying the Rankine- Hugoniot conditions 
at the mean position of the shock wave. 

Carrier (1949a) and Van Dyke (1953) first investigated the response of the inviscid flow 
field to time-periodic oscillations of the wedge. Their interest was in the surface pressure 
distribution and the resultant forces and moments, and their relevance to oscillating airfoils. 
The related problem of freestream turbulence amplification caused through interactions 
with shock waves has been considered by Anyiwo and Bushnell (1982), Hussaini, Collier 
and Bushnell(1986), Meadows, Kumar, and Hussaini (1991), and Kumar, Bushnell and 
Hussaini (1989) who confirmed the important result found in Moore (1951), Ribner (1954) 
and McKenzie and Westphal (1968) that a pure acoustic, vorticity, or entropy wave, upon 
interacting with a plane shock wave, generally generates all three classes of disturbance 
downstream. The interaction of a shear wave with a detonation wave induced by a wedge 
in a supersonic flow, as well as the response of the oblique detonation to oscillations of 
the wedge was considered by Lasseigne and Hussaini (1992). In studying the interaction 
of a shear wave with the detonation, they assumed that a weak steady sinusoidal vorticity 
wave is obliquely convected through an overdriven detonation attached to a wedge. Their 
concern was with the response, measured by the deviation of the detonation position from 
its unperturbed state and by the vorticity and pressure at the detonation, as the degree 
of overdrive is increased. They compared the flow in the presence of the wedge with the 
unobstructed flow field considered in Jackson, Kapila, & Hussaini (1990). In particular, 
for this special type of disturbance (t.e. steady), they found an infinite discrete spectrum 
of disturbances to exist behind the shock. 

The present study also raises other important, broader issues, notably the stability of 
the shock on the wedge. It is generally accepted that the so-called “strong shock” solution 
(Liepmann and Roshko 1957) is unstable (see Levinson 1945, Carrier 1949b, Henderson 
and Atkinson 1976, Rusanov and Sharakshanae 1980, and Salas and Morgan 1982), whilst 
the proof for the stability of the “weak shock” solution has largely been either numerical 


3 


1111 Pi 111 I III! 


(Rusanov and Sharakshanae 1980, and Salas and Morgan 1982) or subject to some restric- 
tions. The work of Henderson and Atkinson (1976) considered just finite length wedges to 
“avoid unbounded velocity downstream” (sic) whilst Carrier (1949b) did not worry about 
convergence in the large” (sic) when considering series solutions. One of the aims of this 
paper is to place the evidence for the stability of flows in which the flow downstream of 
the shock is supersonic on a much firmer footing. Here we shall remove Henderson and 
Atkinsons’s (1970) restriction and pay particular attention to the convergence of series 
solutions. 

The layout of the paper is as follows. In Section 2 we consider the three distinct classes 
of disturbances upstream of the shock, and then in Section 3, the analytic solution for the 
downstream development of these disturbances in the region bounded by the shock and 
the wedge surface is obtained. In Section 4, we examine various aspects of the analytical 
solution. These raise important questions regarding the stability of the shock itself. We 
show that the weak shock solution is stable to these imposed disturbances, provided that 
the downstream flow is supersonic in nature, in line with previous (mostly numerical and 
experimental) studies. We also examine the behaviour of the far-downstream flow which 
suggests an important physical decomposition of the solution. The physical aspects of the 
decomposition is examined in Section 5 and the responses to the various types of imposed 
disturbances are quantified. Our conclusions are given in Section 6 . 


2. Formulation 

Throughout we will denote quantities upstream of the shock by subscript 1 and down- 
stream by subscript 2. We take the wedge surface to make an angle 6 with respect to the 
on-coming flow, with ( x*,y *) the coordinates parallel and perpendicular to the upstream 
base flow (respectively) which has magnitude U * , Mach Number M \ , density p l , and tem- 
perature T* . We assume that the ratio of specific heats 7 is constant. The velocity vector 
is written as Uf(u t v) with respect to (x*,y*) coordinates. The density field is written as 
p\p, pressure as p* x R*T*p , where R* denotes the gas constant and the temperature field 

as T*T. 
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The Rankine-Hugoniot relations provides a link between conditions upstream and 
conditions downstream of the shock, and in particular leads to the following classical 
result in gas dynamics (e.g Liepmann and Roshko 1957), 

tan (/3 - 9) _ U 2 _ Pj_ _ (7 - l)M t 2 sin 2 ft + 2 
tan/3 ui p 2 (7 + l)Afj 2 sin 2 ft 

where ft is the angle between the shock and the wedge centerline, and u\ and U 2 are the 
non-dimensional velocity components perpendicular to the shock. 

It is well known (Liepmann and Roshko 1957) that for a given upstream Mach number 
Mi and 9 < 9 rnax {M\ ) equation (2.1) admits two possible solutions for ft: the larger value 
corresponding to the so called “strong shock” solution, which is characterized by subsonic 
flow behind the shock in all cases and the smaller value corresponding to the so-called 
weak solution which is characterized by the flow being generally supersonic behind the 
shock, except for a small region of subsonic flow for 9 close to 9 = 9 max . For 9 > 9 max , 
no attached-shock solutions to the problem exist. Attached shock solutions all have the 
property of uniform downstream flow directed parallel to the wedge surface. 

Small amplitude disturbances upstream of the shock may be classified into three dis- 
tinct classes (see McKenzie and Westphal 1968, for example). Taking e to be the measure 
of the amplitude, and hence a small parameter, we have: 

(i) Acoustic waves: These are characterized by having a pressure perturbation and the 
resultant perturbations in velocity, density and temperature; however, the waves carry 
no change in entropy and have no vorticity. In this case the disturbance field upstream 
of the shock is given by 


p=l + eE + 0(e 2 ), 

(2.2) 

ot\eE ^ 

7 M?(a.+«) ( °’ 

(2.3) 

7 M?( a.+u,) + ° (€ >’ 

(2.4) 

P= H + 0(e 2 ), 

7 

(2.5) 
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T= 1 — 


e(7 - 1 )E 


( 2 . 6 ) 


7 

where E is the normal mode exponential 


+ 0(e 2 ), 


E = expjiaq X] + ia 2 yi + iut}, 


(2.7) 


and 

w ± + a 2) 1/2 ( 2 - 8 ) 

is the frequency of the disturbance for the given wavenumbers. The coordinates aq 
and ?/i are parallel and perpendicular to the upstream flow respectively, suitably non- 
dimensionalized (this can be accomplished by using one of the wavelengths of the 
disturbance as a characteristic length scale, e.g. setting aj , say, to unity). The choice 
of frequency oj and also the nature of the solution (2.2)-(2.6) arises from solving the 
irrotational flow problem upstream of the shock. We shall refer to modes correspond- 
ing to the negative sign in (2.8) as so-called “fast modes”, and to the positive sign in 
(2.8) as the so-called “slow modes”. 


(*'*) Vorticity waves: These are characterized by having no density, temperature or pressure 
components to the disturbance to 0(e). As such, the upstream field can be written 


together with 


u = 1 + eE + 0(e 2 ), 

(2.9) 

v = -e—E + 0(e 2 ), 
« 2 

(2.10) 

P,T, p= 1 + 0(e 2 ), 

(2.11) 

u = -e*i, 

(2.12) 


with E having the same definition as in (2.7). 


(in') Entropy waves: These are characterized by having no pressure or velocity components 
to the disturbance to 0(e). As such, the upstream field can be written 

p = 1 + eE -(- 0(e 2 ), (2.13) 
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T= l-eE + 0(e 2 ), 
u,p= 1 + 0(e 2 ), 


(2.14) 

(2.15) 

(2.16) 


« = 0(e 2 ), 

together with u> defined as in (2.12) and E defined as in (2.7). 

In the following section we go on to consider the interaction between these waves 
and the shock attached to the wedge. Fortunately, in spite of their distinct features, the 
analysis is fundamentally the same in all three cases, and we shall see how just one of 
these waves upstream of the shock generally produces a combination of all three modes 
of disturbance downstream of the shock; however, unlike the situation where no wedge is 
present, the downstream disturbance is not restricted to a single wavelength. 


3. The Solution Downstream of the Shock 

Downstream of the shock, the flow variables may be written 

P = P2 + ep + 0(e 2 ), (3-1) 

p = p 2 + ep + 0(e 2 ), (3.2) 

T = T 2 +ef + 0(e 2 ). (3.3) 

The non-dimensional coordinates parallel and perpendicular to the wedge are taken as x 2 
and y 2 respectively whilst the velocity components in the x 2 and y 2 direction are written 
as U 2 + eu + 0(e 2 ) and ev + 0(e 2 ) respectively. The governing equations at 0(e) may be 
written 


Pt + U 2 p X2 + P2U X2 + p 2 Vy 2 — 0, 

(3.4) 

p 2 {u t + U 2 u X2 } + ,, 2 p l2 = o, 

7 Mf 

(3.5) 

P2{vt -1- + ^j^2Py2 ~ 

(3.6) 

p2{Tt + U 2 T x j} {pt + U 2 p X2 ) — 0, 

7 

(3.7) 
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p = P2T + Tip. 


(3.8) 


The general solution procedure is based on that followed by Carrier (1949a) and 
Van Dyke (1953). It is found convenient to first split the solution into two components, 
corresponding to acoustic waves and vorticity waves. Specifically, we write 


u = (f) X2 + E y2 , (3.9) 

v = <f>y 2 ~ E l2 , (3.10) 

where <j> — <f>(x 2 ,y 2 ,t) represents the acoustic mode and E — E(x 2 ,y 2 -,t) represents the 
vorticity mode. 

Substituting these forms into (3.5)-(3.6), differentiating and combining appropriately 
yields the following equations 

V J W,+^ 1 + ^} = 0. (3.11) 

\7 2 {E t + U 2 E X2 } = 0. (3.12) 

Solving (3.11) and (3.12) implies the introduction of two arbitrary harmonic functions (see 
the footnote in Hui 1969). Equations (3.11) and (3.12) taken together represent a sixth 
order system which is derived from the second order system (3.5)-(3.6), and therefore 
without loss of generality, we set these two arbitrary harmonic functions to zero and note 
that all physical conditions at the shock and the wedge boundary are satisfied by the 
resulting solution. Using (3.4), (3.7) and (3.8), it is then possible to show that (f> satisfies 
the following equation 

VV = —z{4>ti + 2U 2 (f> X2 t + (3.13) 

at 


with 


where 


p = r(^t + U 2 <f>x), 


(lg — 


T ^ 1 / 2 

±2 

Mi 


(3.14) 


(3.15) 
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is the nondimensional sound speed behind the shock. A further quantity is also required 
to describe the problem fully, namely ip'(y 2 ,t), the displacement of the shock wave. 

It is possible to write the general solution for <f> and E which is bounded at the tip of 
the wedge as 


where 


E = e 


£ 00 

' Yj { a » cosh(i/02) + ^ sinh(^02)} 

(3.16) 

t/=0 


...... »w*2 iufV 2 Xa i 00 

(3.17) 

i/=0 


. . iw 00 

(3.18) 

0 


II 

(3.19) 

a; 

(3.20) 

w 

w = — , 

a s 

(3.21) 

i= \j x 1 - p 2 , 

(3.22) 

r 2 =x 2 ~p 2 y 2 , 

(3.23) 

tanh 82 = , 

X 2 

(3.24) 

A = cot (P — 9). 

(3.25) 


In deriving (3.16)-(3.18), we have assumed that the harmonic time dependency upstream 
of the shock is transmitted unchanged through the shock, an assumption that is likely to be 
accurate provided that |e| << 1 and the wavelength of the disturbance is large compared 
to the shock thickness. The boundedness argument at the wedge tip (r = 0) precludes the 
Y u (z) family of solutions. 

If we impose the impermeability condition on the wedge surface, then the b„ coeffi- 
cients must all be set to zero; if we demand that the shock remains attached to the wedge 
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tip then do = 0, and if we assume that there is no transient behaviour (i.e. only the 
sustained response to the imposed disturbance is present) then «o = Co = 0 also. 

The key useful results are 

iust— v°° > A ( Uo Uo 1 

p = — - e u 2 ~ a ’ y k < 2 ia v -) a v - \ a v +\ \ cosh(i/# 2 )^(kr), (3.26) 

2 a 9 “ L o>s o s J 


i w y 2 A a ^ OO 


U = t 2V 2 


« , iwTJ *) X n OO 

1 lw< ~ r/i-.g 
_ 2 C 


>; oo ^ \ 

aJ + 2^ C,/+1 ~ 2 ^ C>, ~ 1 ) Jl, ^£ y i) 

x 2 °°^ r jj ^ 

^ £ < 2i — a,, + ay_i — a u + 1 > cosh^t^) J„(fcr), 

,,-n l ° s J 


. iwX 2 00 ■ f 32 


(/=0 


t/2 


/ 2 C/ 2 (t/|_a2) \ «3 

1. iwUoJi OO 


+ -e ^ a ' X;M{v.+«rf.}dnh{^m 


«/=0 


(3.27) 


(3.28) 


where a_i = c_i = <f_i = 0. The set of solutions above now involves three sets of coeffi- 
cients, namely {a„},{cy}, and {dy}. Although these solutions satisfy the impermeability 
constraint on 62 = 0, we have yet to impose the unsteady Rankine-Hugoniot conditions 
(McKenzie and Westphal, 1968) on these perturbation terms at the shock, which serve to 
determine the aforementioned set of coefficients. 

The Rankine-Hugoniot conditions on the shock may be written in the following form 


it 2 + u 2 — = (1 — )R' + —{u\P\ + Uj), 

P2 P'2 P2 

, , 1„ h , a l P 2 , , - , Pi 

1/7 d - ~ ^2 d~ _ ■ — Wi ”f“ r .O/\p\ + _ 

2 P2 27^2 P2 2 2'fM^ui 


/ a-l P2 , a 2 p 2 „ u 1niV , ui«', 

u 0 - r rr b 


2 (7- l)u 2 p2 ' (7- 1)«2P2 

Pi + Pi 


= (1 - ^-)R' + ^ 

U 2 U2 


Here we have written 


(7 — l)Mfu 2 (7 — l)Mfv .2 ’ 

U 2 = + (ui - «2)^- 

R' = vi K + tfli 


(3.29) 

(3.30) 

(3.31) 

(3.32) 

(3.33) 
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where a denotes distance along the (unperturbed) shock. We have taken (u', v ') to be the 
velocity perturbations perpendicular and parallel to the shock, respectively, and (u, v ) the 
basic flow perpendicular and parallel to the shock respectively. 

The four conditions (3.29)-(3.33) may be combined to eliminate the downstream per- 
turbation density p' 2 . It is then possible to insert the expressions for p , u , 5, ?/>', all evaluated 
on the undisturbed shock surface, into these resulting equations to yield the following sys- 
tem of equations at each value of v > 0. 


cos(/3 - 6) - ^£c v -i] 

[2«— a v — a „+ 1 + a„_i] cosh(i/0 o )l 
z a s J 


+ sm(ft - 9) | + a„ +1 ]sinh(V0o) + i 

-(Si - a 2 )l-^<4 + i«^+i - <4-i)]| = R?\ 

Z J 

{ b\ } { sin ^ " *) \}~ iX % Cu + ¥ Cl,+1 - 

~ 2 ^* ~ a "+ 1 + Oi'-i] cosh(i/0 o ) j 

-cos(fi - 9) | + Oj/-)-!] sinh(^^o) + *7? 2 ^-c„ j j 

{bI} Sin ^ ~ ^ + ^(d„ +1 - «?„_!>] + ift 2 a s d v 

M 3 \/l ro . E7 2 , 

M 2 -f* ^ a„_j ~ cosh(t'Po) > 


+ 

+ 


{ 


i? 


_ i 

» fli 3) 


Here we have written 


tanh So = fl tan(/? — 0), 


X ’ = 2 - 


* = 


2^2 ’ 


(3.34) 


(3.35a, 6) 


(3.36) 

(3.37) 

(3.38) 

(3.39) 
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(3.40) 


B x = 1 + 


a. 


(7-I) 


~o ’ 
ui 


- si ) (1 ■ 


7 a * 


(3.41) 

ft = , - ■ (3 42) 

(7 - 1 )u 2 

The determination of ,R?\ and is routine once the upstream disturbance is 

prescribed, but for brevity these are omitted. These quantities arise from the forcing caused 
by the upstream disturbance which is proportional to the exponential (2.7) evaluated at 
the shock position. This exponential is expanded in terms of the Bessel function by using 


_ * y * „»sin ek£y2 

1 — 

> “ 

1 + £|/,0 

(cvi cot ft + a 2 )(cos 8 + A sin 8) + 


where 


8 = sin 1 


exp[i(oq cot ft + a2)( cos $ + A sin 8 ) 1 / 2 ] — e 1,7 “• e 

• „ 00 —ii/Q _i f 1 \v p — ivQ 

J 2 - — MHm), 


. 00 t i/O 1 f 

-*7 Tt77 yi e +1- 1 ) e 

= e 2 • 

i /=0 


(3.43) 




(3.44) 


and 6„ )0 is the Kronecker delta function. Then a u + 1, 6„+i, and c„+ 1 are determined from 
equations (3.34), (3.35) which are essentially recursive in form. 


4. Analysis of the Series Solution and Far-Downstream Behaviour 

In Figure la, we show %{pe~ iu,t } (solid line) and 3{pe -,u,t } (broken line) both evalu- 
ated along the wedge surface for the particular case M\ — 5, 8 = 25°, oq = 1, and £*2 — 0) 
taking the positive sign in (2.8) (the so called “slow mode”). The corresponding shock po- 
sition %{xi>'e~ iwt } (solid line) and Q{^e- iut } (broken line) with distance along the shock, 
7. 2 = 1/2(1 + A 2 ) 1 / 2 is shown in Figure lb. These results indicate a complex structure to 
the solution; both quantities show the expected oscillatory behaviour, although the wedge- 
surface pressure shows indications of amplitude decay as r 2 -» 00. Shown in Figure 2 is the 
response to the so-called “fast mode”, i.e. the negative sign in equation (2.8) is used. The 
wedge-surface pressure is shown in Figure 2a whilst the shock position is shown in Figure 
2b; as before, real components are indicated by solid lines and imaginary components by 
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broken lines. The expected oscillatory behaviour is shown, but in this case there is little 
evidence of amplitude modulation of the wedge-surface pressure. These two examples indi- 
cate the importance of ascertaining the correct downstream behaviour of the solution given 
by the infinite sums (3.26-3.28). Unfortunately, it was found that as i % 2 (or £ 2 ) became 
large, the computations for xp 1 (and p) become extremely sensitive to the number of terms 
of the series taken (although we stress that all results presented are carefully checked for 
accuracy, in particular for insensitivity to truncation of the number of terms in the Bessel 
series (3.16)-(3.18)). Cursory inspection of the coefficients a u , b„, and d w indicated that 
the numerical difficulties described result from these coefficients often becoming extremely 
large as u — ► 00 , leading to numerical difficulties associated with the effects of roundoff. 

The same general trend of oscillatory, complex behaviour seen in Figures 1 and 2 was 
also seen in the results of Lasseigne and Hussaini (1992) relating to detonation waves and 
moving wedges. However, there are two, not unrelated, issues that remain unanswered and 
affect the physical interpretation of the results. These two issues are: (i) the nature of the 
coefficients a„, c„ and d v as u —* 00 and (ii) the precise nature of the flow far downstream 
of the shock. We look at these aspects in some detail and discover a number of interesting 
subtleties. 

Superficial inspection of the coefficients indicates general growth as u — * 00 . Indeed, 
this was observed even when the inhomogeneous forcing terms R\}\ R^\ and R[?\ were 
“turned off’ after v = 1, in order to mimic a homogeneous system (although in this case 
it was found that d„ — > 0 as v —* 00 ). We now therefore consider the homogeneous system 
for (3.34), (3.35 a,b) in the limit as u — * 00 . This leads to a system of difference equations, 


and we therefore surmise that the coefficients take on the following form as u — » ex? 

c„ « C 0 K\ (4.1) 

d v « D 0 K\ (4.2) 

whilst because of the nature of the occurrence of the a„, we must have 

a„ « A Q K v t~ v9 \ (4.3) 


where K is an eigenvalue whose precise value is extremely important since it is this value 
that would be expected to determine the ultimate downstream behaviour of the disturbance 
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field. If | A' | > 1, then using the generating function for Bessel functions yields the following 
result 

OO OG 

Y c v Jv{z) » Y, C*K v J u {z) « C 0 e^ K -'' K)z , \z\ - oo. (4.4) 

V— 1 t/=l 

This then predicts exponential growth as |z| — ► oo if 9?{z(A' — 1/K )} > 0. Therefore, we 
proceed to the actual determination of K. 

Substituting (4.1)-(4.3) into (3.34), (3.35a,b), and discarding inhomogeneous terms, 
leads to the following three equations to determine K : 


cos(/3 - 9) + i (K - ijif-'jCo 

— 7 [2i — - e~ e °K + e e ° AT _1 ]>lo 1 
4 a s J 

+ sin(/3 - 9) | ^p[e- 9 °K + e 6 ° K~ l ]A Q 

+ifi 2 y7-Co - (B, - fi 2 )[-— + h'A' - ifjr-'lZJo} = 0, 

U 2 2» 2> J 


(4.5) 


- y[2t— - Ke~ 9 ° + K~ l e e °)A Q 
4 a s 

-co S (,9-9){i9|A> (46a>J) 

+ j ^ } sin(/J - 9) - A’ 1 )) + i/3 2 a,]Do} 

= 0. 

This system may conveniently be expressed in the form of a sixth-order linear gener- 
alized eigenvalue problem by writing 


Ao = KAq , 


Co = KCo , 


(4.7) 


Do = KDo, 
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and solving for the eigenvector {A 0 ,A 0 ,Co,Co,D 0y Do}. This was accomplished numeri- 
cally and yielded the following general trends: (i) Of the six roots for K of the eigenvalue 
problem, four were imaginary, and the remaining two were complex twins K = A'o and 
K — —c.c.{Ko} (indeed, the nature of these roots can be confirmed by inspection of the 
system (4.5), (4.6a,b); (ii) the imaginary roots all had 0 < | A'l < 1; (iii) the complex pair of 
roots had \K\ > 1. This last point strongly suggests the flow is structurally unstable since 
there exists homogeneous solutions for <j> and E which grow exponentially as the distance 
from the wedge tip increases; however, further detailed analysis given later shows this to 
be untrue. 

Returning now to consider the nature of the imaginary roots, Figure 3 shows the 
variation of these roots with wedge angle 9 at M\ = 1.5, 3, and 5. In addition to the 
trends outlined above, the roots have the features that K — ► 0 as 9 — > 0, whilst as 9 — * 9 m , 
IA'1 — > 1, where 9 m is the limit of the envelope of supersonic flow downstream of the shock. 
This then, in some way connects with the results of previous research on the stability of 
the strong shock solution. We now consider the nature of the complex roots. Fortunately, 
it is possible to obtain precise analytic information on these roots aided by the observation 
that the eigenvectors from the computation reveal that Do = Do = 0 (to within machine 
precision). If we use this result, then the groupings shown in (4.5), (4.6) would suggest that 
an eigensolution is possible if the following three equations are simultaneously satisfied: 


+ Uk - Uk-^Co - j[2i— - e~ d ° K + e e ° K~ l ]Ao = 0, 
U 2 2 2 4 a 3 


][e~ 8 °K + + ip^Co = 0 , ( 4 . 9 ) 

4 U 2 

2 i - —e~ 6o K + — e e ° K~ l = 0 . ( 4 . 10 ) 

Cl ff Clg 

We find this is indeed possible and the solution is 

K= e ^-{ia,±[Vl-ci\)'l 1 ), ( 4 . 11 ) 

u 2 

Co = ±iAo/2. ( 4 . 12 ) 
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The value of K, given by (4.11) has \K\ — e 9 ° >1, suggesting that the flow is always struc- 
turally unstable in the sense that disturbances will always grow downstream. However, 
it turns out that in fact this homogeneous solution represents a structurally stable flow. 
This stability is demonstrated by defining the following functions which arise naturally as 
a result of the analysis above 


E = C 0 e 




U 2 (Ul 




(4.13) 


together with 


iut- 


Uixi oo 


4> = A 0 e V S-: Y, {Ke 9 ^y J„(*r); (4.14) 

v~ — OO 

E and <f> are both harmonic and related through the Cauchy-Riemann equations, with 


4>X 2 — Ey 2 , <fiy 2 — E X2 . 


(4.15) 


As a direct consequence of (4.15), the exponential growth of the two components to this 
homogeneous solution cancel each other when computing physical quantities such as the 
pressure p or the velocity (u,5), and therefore, the flow calculated from the summations 
(3.26)-(3.28) is bounded for all values of ( 22 , 2 / 2 )- 

Although, the growing homogeneous solution for the coefficients does not lead to a 
physically unbounded solution, due to exact cancellations when calculating physical quan- 
tities, their presence greatly complicates the computation of the summations (3.26)-(3.28) 
for calculating the flows since the relevant physical quantities result from the remainder of 
summing three or more, numerically very large quantities, to yield a much smaller quantity. 
The numerical difficulties in obtaining accurate solutions associated with this approach are 
obvious. We note (again) that all results presented have been carefully checked, to account 
for this phenomenon. 

The final question that arises is the nature of the downstream response to sustained 
acoustic waves, entropy waves, or vortic.ity waves. We may expect that the solution down- 
stream will remain oscillatory, with no amplitude decay, although this is not always the 
case. The different downstream behaviours of the acoustic disturbances can be understood 
by firstly examining the solution in the absence of the wedge (see McKenzie & Westphal 
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1968). Since the frequency of the disturbance downstream of the shock is fixed by the fre- 
quency of the disturbance upstream of the shock and the wavenumber of the disturbance 
tangential to the shock is also fixed by the deflection of the shock produced by the upstream 
disturbance, there is a finite range of wavenumbers tangential to the shock that exclude the 
existence of plane acoustic waves which propagate downstream. If the imposed tangential 
wavenumber lies within this range, then the acoustic disturbance behind the shock (in the 
absence of the wedge) exhibits exponential decay in a direction normal to the shock. If the 
tangential wavenumber lies outside of this range, then a plane wave propagates at an angle 
6 P to the shock. This angle and the range of excluded tangential wavenumbers depend 
upon the relative orientation of the shock and the upstream base-flow and disturbance 
conditions. 

The presence of the wedge changes the nature of these solutions drastically. In the 
case of the tangential wavenumber being within the excluded range, an acoustic field is 
generated that decays algebraically (rather than exponentially) owing to the requirement 
that the shock remain attached to the apex and the requirement that disturbances have 
zero normal velocity at the wedge surface. If the tangential wavenumber is outside of 
the excluded range, then the solution in the presence of the wedge exhibits two types of 
behaviour, dependent upon the angle of propagation of an acoustic disturbance in the 
absence of the wedge. If the angle of propagation 9 p is greater than the angle between the 
shock and the wedge /? — 6, then the pressure disturbances generated at the shock would 
intersect with and reflect from the wedge surface. Thus, there would be a non-decaying 
pressure field everywhere between the shock and the wedge surface with a superimposed 
algebraically decaying pressure field owing to the requirement of an attached shock and 
zero velocity normal to the wedge surface. If the angle of propagation 9 P is less than the 
angle between the shock and the wedge surface, then the pressure field is divided into two 
regions by a ray emanating from the apex and parallel to the direction of propagation 
of the pressure disturbance in the absence of the wedge. Between the shock and this ray 
there is a non-decaying pressure field, with the aforementioned superimposed non-decaying 
component. Between the ray and the wedge surface the pressure field is algebraically 
decaying and there is a component with sustained oscillations. This is in contrast to the 
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solution in the absence of the wedge and is explained by the presence of the wedge that 
serves to terminate the shock (which is assumed of infinite extent in both directions in the 
absence of the wedge), and there are no pressure disturbances generated by the shock to 
fill this region. We will refer to this region as the shadow region. 

The above mentioned downstream limits can be mathematically understood by again 
considering the behaviour of the coefficients a„, c„, and d*, as v — * oo when the inho- 
mogeneous terms R/J , i?t 2 \and are included. For this, the nature of 8 in (3.44) is 
important. If 8 is real, then the wavenumber of the imposed disturbances lies within the 
excluded range derived from the theory of McKenzie and Westphal (1968), the coefficients 
within the summation (3.43) neither grow nor decay, and the combination of coefficients 
within the bracket in equation (3.26) all decay at a rate proportional to e - "^ 0 . Thus when 
reconstructing the pressure along any ray $•_, = constant , algebraic decay is found since 
# 2 must lie on the interval (0, # 0 ) with zero corresponding to the wedge surface and #0 
corresponding to the shock. For 8 2 = #0 the bracketed combination of coefficients does 
not decay leading to the appropriate oscillatory behaviour necessary to match conditions 
at the shock. If 8 is complex, then the condition 

(oq cot (3 + O 2 )( cos 8 + A sin 8) + 

M 

is satisfied. We also note that for equation (3.43) to be satisfied then sin# must remain 
real, which requires that 

8 = ±7r/2 — i8i, (4.18) 


> 1, (4-17) 


where 


(oq cot /3 + G; 2 ) ( cos 8 + A sin 8) + 


Notice 8i is related to the angle between the direction of the shock and the direction of the 
acoustic disturbance 8 p through tanh 8, = tan 8 p . For imaginary values of 8 the coeffi- 
cients of (3.43) behave as (±ie 9i ) l/ , as u — > 00 , which increase exponentially in magnitude 
as do the corresponding coefficients R.[}\ R^\ and R^\ If #; is greater than 8q, then the 
generated or transmitted acoustic waves intersect the wedge surface. The physical conse- 
quence of this intersection is seen by considering the bracketed combination of coefficients 
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in (3.26) along with the factor cosh(n#2)- The product grows exponentially and alternates 
between real and imaginary. The relation (4.4) shows that there is a sustained oscillatory 
solution of the pressure for all and j/2- For 0 < 6i < 6o, the aforementioned product 
of coefficients is exponentially growing when 62 > Oi which shows that there is a region of 
sustained oscillatory behaviour of the pressure field near the shock, and the aforementioned 
product of coefficients is exponentially decaying for 6 2 < 6i which shows that there is a 
region of algebraic decay of the pressure field near the wedge surface and that there are 
no sustained oscillations near the surface. This overall behaviour of the pressure field has 
been verified by direct computation of the summation (3.26) along rays 62 — constant for 
various cases. 

With the previous analysis in mind and restating that one goal of this analysis is to 
discover how these solutions might affect the boundary layer receptivity problem, we focus 
attention on the pressure evaluated along the wedge which generally may naturally be 
decomposed into four parts 


e -i*t p = Poe -'^f+ikcosm-6 0 )x 2 

+ P,x-^ 2 e-^-^ 

+ P 2 X2 l / 2 e~ ( ' ^ 2 +“* + Q(x 2 ). 


(4.20) 


The first term on the right-hand-side is non-zero only when 0i > 60 and is evaluated 
numerically by considering the large u behaviour the coefficients |o„}. This term is con- 
verted to an infinite sum of Bessel functions and then subtracted from the left-hand-side 
where the pressure has the representation given by the summation (3.26). The Fourier 
Sine Transform of the resulting difference is evaluated and the strengths of the inverse 
square root singularities provide the amplitudes Pi and P2 for the second and third terms. 
Physically, the first term represents the transmitted acoustic wave (or generated acoustic 
wave if the imposed disturbance is an entropy or a vorticity wave) when it intersects the 
wedge surface, and the second and third term are the far-downstream behaviour of slow 
and fast moving acoustic waves associated with reflections off the wedge surface and regen- 
eration of new acoustic waves at the shock, as well as new acoustic waves generated by the 
requirement that the shock remains attached to the wedge tip. The first three terms are 
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1/2 

the leading terms in the limit x 2 —* oo, and the fourth term decays faster than x 2 ' and 
is considered to be the local field. The Fourier spectrum of the fourth term is important 
since it represents the generation of all wavelengths by the response of the attached shock 
to a single wavelength input; the spectrum of the local field is bounded for all wavenumbers 
since the local field decays relatively quickly as 12 00 . 

If0, < 0 O then the far- downstream behaviour of the pressure along the surface changes. 
The wavenumber of the imposed disturbance is outside of the allowable wavenumbers for 
propagating plane waves or the generated or transmitted plane- wave pressure disturbance 
does not intersect the wedge, and it is found that Po is zero. The amplitudes P\ and P 2 
are non-zero and therefore the pressure has algebraic decay far downstream. 

The vorticity, found by computing O = V 2 £, can be similarly decomposed into four 
parts if 9i > 0 

__ iw a \ \ y 2 

e“ ,u ' <+ ^ i O = c _ u 2 (i/|-.J)^ (je .'fcco»h# i »2 

+ fi 2 y 2 " I/ 2 e i ^ +, T}+«(y 2), 

where we have assumed that the calculations take place far enough downstream in x 2 such 
that the behaviour represented by (4.21) is an accurate expansion valid for large y 2 . The 
first term represents a generated or transmitted vorticity wave and the next two terms 
represent additional vorticity generated by the moving shock as a result of requiring that 
the shock remain attached to the wedge apex and zero normal velocity on the surface. That 
is, the algebraically decaying component of the pressure field generated by the presence of 
the wedge interacts with the shock wave through the Rankine-Hugoniot relations which 
indicate that the shock must be distorted by this component of the pressure field, and 
hence algebraically decaying vorticity modes represented by and fl 2 are generated 
by the curvature of the shock. A numerical subtlety exists when trying to calculate the 
decomposition (4.21) when 9 t > 0. The coefficient of the vorticity depends on cancellations 
involving five very large terms (for the pressure decomposition there were only three very 
large terms); thus, it is not possible to accurately determine for some cases. If 6 t = 0, 
then it appears to be not possible to calculate a decomposition of the vorticity into the form 
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given in (4.21). The coefficients of the vorticity expansion neither grow nor decay making 
the numerical determination of flo difficult if not impossible. In the absence of the wedge, 
a single wavelength vorticity mode would still exist for this case, and the pressure field 
would be exponentially decaying. With the wedge present, there exists a small algebraic 
component of the pressure field. 

The results relating to the decompositions given in (4.20) and (4.21) are given in the 
next section. 


5. Results 

The response of the flow downstream of the shock is generally complicated, although 
the analysis of the previous section does enable us to quantify the response, at least in 
the far downstream limit. To this end, we find it useful to consider the sine-transform of 
Q(x 2 ), that is 

pOO 

Q*~ Q(x2)sin(fiX2)dx 2 . (5.1) 

Jo 

Fortunately, this quadrature can he carried out semi-analytically using the following sine- 
transform 


fOO 

/ sin(iJ,X2)e tqX2 J l/ (kx 2 )dx 2 = 

Jo 


, | \A 2 -(p+s) 2 +*(/*+«)) 

A " -(#*+ «> 2 

j | - 0 * - «) 2 - >0 - ?)) 


(5.2) 


2ik v 


\ A 2 - - «) 2 


The inverse square-root singularities lead to the second and third terms of (4.20) and 
the evaluation of (5.1) involves a summation of terms derived from (5.2) with coefficients 
determined numerically. The results of this process are shown graphically in Figures 4-8 
where 3?{Q*} is denoted by solid lines, whilst £l{Q*} is denoted by broken lines. 

The spectrum Q* for the input conditions M\ = 5, <*1 = 1, <*2 = 0 and 0 = 25° 
(slow mode) is shown in Figure 4. For this set of input conditions, inequality (4.17) is not 
satisfied and therefore Pq in (4.20) is zero. This spectrum is bounded for all wavenumbers 
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but has square-root behaviour at /j = , the wavenumber locations of the slow and 

fast acoustic disturbances behind the shock. This square-root behaviour remains after 
subtracting the inverse square-root singularity to form the second and third terms of (4.20). 
For wavenumbers between these two points, a plane wave with the fixed time frequency u 
would have to travel faster than the sound speed relative to the mean flow, and so these 
represent pressure disturbances that decay exponentially in x 2 . The local field is comprised 
of a Fourier integral over both classes of disturbance. It should be noted that the large 
wavenumber disturbances are “slow modes” whilst the small wavenumber disturbances are 
“fast modes” . The importance of Q* is that it demonstrates that far downstream the energy 
density of the acoustic response is spread over all wavenumbers and is not concentrated 
at an isolated wavenumber as is the energy density of the acoustic mode in the absence of 
the wedge. In particular, the energy density is significant at large wavenumbers, i.e small 
wavelengths, and these disturbances generated by the shock-body interaction might lead 
to receptivity as they interact with a growing boundary layer on the wedge surface. 

Figure 5 shows the spectrum for the fast mode with the same input conditions as in 
Figure 4. For the fast mode, inequality (4.17) is satisfied so that P 0 is non-zero. To indicate 
the presence of this transmitted plane wave, a circle has been added to Figure 5 at the 
wavenumber location of the transmitted wave. It should be noted that this wave number 
location is not quite at the same wavenumber as the singularity of the Fourier transform 
which leads to the second and third terms of the expansion. Figure 4 and Figure 5 appear 
to be quantitatively similar, although a shift to smaller wavenumbers is forced by the 
changed frequency of the imposed disturbance, a result of choosing the opposite sign in 
equation (2.8). 

In Figure 6, the spectrum of a fast mode with larger incoming Mach number is shown. 
The input conditions are M\ = 10, ac\ = 1, Q.% = 0 and 9 = 25° (slow mode). This 
set of input conditions also falls within the envelope of flows which lead to sustained 
downstream acoustic waves. It is seen that the energy density of the local field has a 
strong Mach number dependence, as do the amplitudes of the first three terms of (4.20) 
which will be explored in Figures 9-12. Figures 7 and 8 show the spectra for the local 
pressure field Q * , in response to an incoming entropy disturbance (2.13-2.16) and to an 
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incoming vorticity disturbance (2.9-2.12). The input conditions are the same as in Figure 
4 (except for Figure 8, for which we have taken oq = 0 and c *2 = 1)- There is a qualitative 
difference between the local field response to acoustic disturbances and the local field 
response to these non-acoustic inputs. The spectra in Figures 7 and 8 are considerably more 
peaked near the fast-mode singularity with the response near the slow-mode singularity not 
nearly as strong. For the response to both the fast and slow incoming acoustic modes, the 
spectra near both singularities were broad with the strength of the slow-mode singularity as 
strong or stronger than the fast-mode singularity. Perhaps this is a result of the incoming 
entropy disturbance having zero pressure and velocity components, the incoming vorticity 
wave having zero pressure and density components, whilst all components of the incoming 
acoustic disturbances are non-zero, thus leading to a much more complicated interaction 
with the attached shock (especially near the tip of the wedge). 

The amplitude of the transmitted or generated (sustained) acoustic disturbance P 0 , 
the amplitude of the algebraically decaying slow mode P \ , and the amplitude of the alge- 
braically decaying fast mode P 2 are shown graphically in Figures 9-12. It can be shown 
that Pa and P 2 are real quantities whilst P\ is strictly imaginary. 5? { P 0 } is denoted as a 
solid line, ^{Pj} is denoted as a line with small dashes and 3 ?{P 2 } is denoted as a line 
with large dashes. 

Figure 9a shows the response for the an impinging slow acoustic mode. Over a large 
range of Mach numbers of the incoming flow M \ , there exists no downstream sustained 
acoustic waves, i.e Pq = 0; however, both of the algebraically decaying modes are non-zero, 
and the amplitudes increase considerably as the Mach number increases. As mentioned in 
the introduction, this is significant since in the absence of the wedge the theory (as verified 
by numerical computations) predicts that the acoustic response will be exponentially de- 
caying, but by requiring that the shock remains attached to the wedge tip, it is seen that 
the response decays like the inverse square root (together with a local field which also does 
not decay exponentially). Figure 9b shows the details of the transition region (in Mach 
number) between the decaying pressure response and the sustained pressure response. The 
amplitudes of Pj and P 2 jump dramatically when Po is non-zero. Physically, this implies 
that there are two parts to the algebraically decaying modes. The first part is a result of 
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the requirement that the shock remain attached as seen when Po — 0; the second part is 
a direct result of the sustained acoustic disturbance reflecting off the wedge surface and 
interacting with the shock to generate more acoustic disturbances. As M\ decreases, the 
critical wedge angle (above which there is no mean flow solution) decreases until it reaches 
the fixed wedge angle of Figures 9-12 (i.e. 6 = 25°). The amplitude of the transmitted 
wave remains finite in this limit; however Pi and P 2 increase sharply in magnitude, this im- 
plies that the solution is becoming increasingly unstable as the critical angle is approached 
(as would be expected). In Figures 10-12, the response to an incoming fast acoustic wave, 
to an incoming entropy wave and to an incoming vorticity wave, respectively, are shown; 
in all cases 6 = 25°, Mj = 5, <*i = 1, and a 2 = 0 except in the case of the vorticity mode 
(Figure 12j where ol 2 ^ 1 and ou — — 0. There is a sustained acoustic disturbance for these 
three cases, and correspondingly, the amplitudes of Pi and P 2 are relatively large since as 
mentioned in the discussion of Figure 9 these two modes have components related to the 
reflection of the acoustic waves from the wedge surface. All three figures demonstrate the 
same behaviour as Figure 9 when the critical angle is approached - the amplitude of the 
sustained mode remains finite whilst the other modes grow. For the responses to the im- 
pinging acoustic disturbance (Figure 10) and to the impinging entropy disturbance (Figure 
11), the algebraically decaying fast-mode is considerably larger than the slow-mode with 
the relative importance of the slow-mode being greater for the response to the entropy 
wave. Comparing the spectra shown in Figures 5 and Figure 7, the energy of the slow 
mode is seen to be in the local fields with relatively more energy for the slow-mode in the 
response to the impinging acoustic wave. In the response to the incoming vorticity distur- 
bance (Figure 12), the amplitude of the slow mode P\ is much larger than the amplitude 
of the fast-mode P 2 . Comparing with Figure 8, this is consistent with the relative absence 
of energy in the slow-mode of the local field. 

For completeness, Figure 13 and Figure 14 show the result of a decomposition of the 
vorticity into the form given by (4.21). Figure 13 represents a response to an incoming 
fast acoustic wave with the same conditions as Figure 10, but with a wedge angle of only 
20°. The change of wedge angle is due to the numerical difficulties of calculating Qo al- 
ready mentioned. Figure 14 shows the response to an incoming entropy wave. General 
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conclusions about the decomposition are difficult to obtain (again related to the numerical 
difficulties); however, we do note the relative insensitivity to the Mach number of the coef- 
ficients except near the critical angle where one of the algebraically decaying components 
becomes quite large. Again, this probably relates directly to the stability of the entire 
solution as the critical angle is approached. 


6. Conclusions 

We have considered the interaction between the shock wave attached to the wedge and 
freestream disturbances. In order to quantify the results, the freestream disturbances were 
taken to be fixed-frequency, single-wavelength, plane-wave disturbances of acoustic type, 
the vorticity type, or entropy type. The behaviour of the flow field downstream of the shock 
in the shock- wedge combination was compared with the behaviour of the downstream flow 
in the absence of the wedge. After rationalizing a number of mathematical and computa- 
tional subtleties, we were able to show precisely how the presence of the wedge affects the 
interaction between a shock and a free stream disturbance. We specifically focused on the 
pressure disturbance along the wedge surface as the analysis of these disturbances directly 
relates to boundary-layer receptivity in the shock-wedge configuration and we also studied 
the vorticity generated by the presence of the wedge as this analysis relates to enhanced 
mixing by the vorticity production of a moving shock. 

The most significant differences in the interaction produced by the presence of the 
wedge are the presence of an acoustic and vorticity field which decays algebraically in 
addition to sustained oscillations produced by the incoming disturbance. It has been 
shown that this field is comprised of two parts: one produced by the requirement that the 
shock remain attached to the apex of the wedge and the other produced by reflections of 
acoustic waves from the wedge surface. 
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